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Abstract. We present radio observations at three fre- 
quencies and contemporaneous optical monitoring of the 
peculiar BL Lac object AO 0235-f 164. During a three- week 
campaign with the VLA we observed intraday variabil- 
ity in this source and found a distinct peak which can 
be identified throughout the radio frequencies and ten- 
tatively connected to the R-band variations. This event 
is characterized by unusual properties: its strength in- 
creases, and its duration decreases with wavelength, and 
it peaks earlier at 20 cm than at 3.6 and 6 cm. We dis- 
cuss several generic models (a "standard" shock-in-jet 
model, a precessing beam, free-free-absorption in a fore- 
ground screen, interstellar scattering, and gravitational 
microlensing) , and explore whether they can account for 
our observations. Most attempts at explaining the data on 
0235-f 164 require an extremely small source size, which 
can be reconciled with the 10^^ K inverse Compton limit 
only when the Doppler factor of the bulk flow is of order 
100. However, none of the models is completely satisfac- 
tory, and we suggest that the observed variability is due 
to a superposition of intrinsic and propagation effects. 

Key words: Galaxies: active - BLLacertae objects: indi- 
vidual: AO 0235+164 - Radio continuum: galaxies 



1. Introduction 

The radio source AO 0235-1-164 was identified by Spinrad 



weeks (e.g Chu et al. |I996| , O'Del l et al. |1988| , Terasr anta 
et al. |1992|, Schramm et al. |1994|, Webb et al. |1988|, this 



& Smith (1975) as a BLLac object due to its almost fea- 
tureless optical spectrum at the time of the observation, 
and due to its pronounced variability. Long-term flux den- 
sity monitoring in the radio and optical regimes have re- 
vealed strong variations and repeated outbursts with large 
amplitudes and timescales ranging from years down to 



paper. Fig. 0). Furthermore, intraday variability in the ra- 



dio (Quirrenbach et al. 1992, Romero et al. 1997), in the 
IR (Takalo et al. 1992), and in the optical regime (Heidt & 
Wagner |l996 , Rabbette et al. 1996) has also been observed 
in this object. In the high energy regime, 0235-1-164 was 
detected w ith E GRET on board of the CGRO (v. Mon- 
tigny et al. 1995), showing variability between the individ- 
ual observations. Madejski et al. ( 1996D report variability 
by a factor of 2 in the soft X-rays during a ROSAT PSPC 
observation in 1993. VLBI observations (e.g. Shen et al. 



1997, Chu et al. 1996, Baath 1984, Jones et al. 1984) re- 



veal a very compact structure and superluminal motion 
with extremely high apparent velocities (perhaps up to 

/3app^30). 

Three distinct redshifts have been measured towards 



0235-1-164 (e.g. Cohen et al. 1987). Whereas the emission 
lines at z = 0.940 have been attributed to the object 
itself, two additional systems are present in absorption 
{z — 0.851) and in emission and absorption {z = 0.524). 
Smith et al. (1977) observed a faint object located about 
2" south of 0235-1-164, and measured narrow emission lines 
at a redshift of z = 0.524. Continued studies on the field of 
0235-1-164 have revealed a number of faint galaxies, mostly 
at a redshift of z = 0.524, including an object located 1'.'3 
to the east and 0'.'5 to the south (e.g. Stickel et al. 1988, 
Yanny et al. 1983). Recently, Nilsson et al. (1996) inves- 
tigated 0235-1-164 during a faint state and found promi- 
nent hydrogen lines at the object redshift oi z = 0.940. 
They note that 0235-1-164 - at least when in a faint state 
- shows the spectral characteristics of an HPQ. Further- 
more, through HST observations of 0235-f 164 and its sur- 



rounding field, Burbidge et al. (1996) discovered about 30 
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faint objects around 0235-1-164 and broad QSO-absorption 
lines in the southern companion, indicating that the lat- 
ter is an AGN-type object. Due to the presence of several 
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foreground objects, gravitational microlensing might play 
a role in the characteristics of the variability in 0235+164, 



as was suggested by Abraham et al. (1993). 

In this paper, we further investigate the radio vari- 
ability of 0235+164, and attempt to determine the most 
likely physical mechanisms behind the observed flux den- 
sity variations. The plan of the paper is as follows. In Sec- 
tion 2 we describe the observations and the data reduction; 
subsequently we analyze the lightcurves and point out 
some of their special properties. In Section 3 we explore 
different scenarios which could explain the variability: we 
discuss relativistic shocks, a precessing beam model, free- 
free-absorption, interstellar scattering, and gravitational 
microlensing. Finally, in Section 4 we conclude with a sum- 
mary of the failures and successes of these models. 

Throughout the paper we assume a cosmological inter- 
pretation of the redshift, we use Hq = 100/ikm/(sMpc) 
and go = 0.5, which gives for the redshift of 0235+164 a 
luminosity distance of 3280 h~^ Mpc; 1 mas corresponds to 
4.2 h~^ pc. The radio spectral index is defined by S,y oc i^". 

2. Observations and data reduction 

2.1. Radio observations 

From Oct 2 to Oct 23, 1992, we observed 0235+164 with 
a five-antenna subarray of the VLAR during and after 
a reconfiguration of the array from D to A. The aim of 
these observations was to search for short-timescale vari- 
ations in several sources. The complete data set will be 
presented elsewhere (Kraus et al., in preparation). In par- 
allel, optical observations were performed in the R-band 



(see Section 2.2). Data for 0235+164 were taken at 1.49, 
4.86, and 8.44 GHz (A = 20, 6, 3.6 cm) every two hours 
around transit, i.e., six times per day. These three sets of 
receivers have the lowest system temperatures and high- 
est aperture effic iencie s of those available at the VLA (see 
Crane & Napier 1989); and data in these bands are less 
susceptible to problems with poor tropospheric phase sta- 
bility than those at higher frequencies. In addition, intra- 
day variability of compact flat-spectrum radio sources did 
appear most markedly in this frequency range in previous 
observations (Quirrenbach et al. 1992| ). During the first 
week, the antennae included in our subarray were changed 
repeatedly due to the ongoing reconfiguration; however, an 
attempt was made to maintain an approximately constant 
set of baselines. Since 0235+164 and the used calibrator 
sources used are extremely compact (cf. VLA calibrator 
list) , the effect of the ongoing reconfiguration on the mea- 
surements is negligible. 

After correlation and elimination of erroneous data 
intervals, we performed phase calibration first. Subse- 
quently, a (one-day) mean amplitude gain factor was de- 



^ The Very Large Array, New Mexico, is operated by As- 
sociated Universities, Inc., under contract with the National 
Science Foundation. 



rived using non-variable sources such as 1311+678, which 
have been linked to an absolute flux density scale (Baars 
et al. 1977, Ott et al. 1994) by frequent observations of 
3C286 and 3C48. After a second pass of editing spurious 
sections of the data, the visibilities of each scan were inco- 
herently averaged over time, baselines, polarization, and 
IFs. Because of the point-like structure of the sources, the 
mean source visibility is proportional to the flux density. 
Eventually, systematic elevation and time-dependent ef- 
fects in the lightcurves were removed, using polynomial 
corrections derived from observations of the calibrator 
sources 0836+710 and 1311+678. 

The errors are composed of the statistical errors from 
the averaging and a contribution from the residual fluctu- 
ations of the non-variable sources 3C286, 1311+678 and 
0836+710. The level of these fluctuations was estimated 
from a running standard deviation of the calibrator mea- 
surements over a two-day period. Over the full three-week 
period, the standard deviations were found to be 0.5, 0.5, 
and 0.7% of the mean value at 1.5, 4.9, and 8.4 GHz re- 
spectively (with no significant difference between the three 
non- variable sources). 

The resulting lightcurves for the three frequencies 
are displayed in the top panels of Fig. H. The mean 
flux densities are 1.57 Jy, 4.05 Jy, and 5.22 Jy for i^ = 
1.49, 4.86, 8.44 GHz, respectively. Therefore, 0235+164 
had a highly inverted spectrum at the time of the observa- 
tions with spectral indices a43GHz ~ 0-80 and ag^QHz ~ 
0.46. 

2.2. Optical observations 

The radio data were supplemented by observations at 
650 nm (R-band filters) taken at the following telescopes: 
0.7 m Telescope, Landessternwarte Heidelberg, Germany; 
1.2 m Telescope, Observatoire de Haute Provence, France; 
1.2 m Telescope, Calar Alto, Spain; 2.1m Telescope, 
Cananea, Mexico. 

Owing to limited observing time per source and 
weather limitations, the optical data are sampled more 
sparsely. They cover the first week of the radio observa- 
tions, leave a gap for ten days, and continue for a total 
of thirty days, thus ending ten days after the radio mon- 
itoring. After the usual CCD reduction process (see e.g. 
Heidt & Wagner 1996), we performed relative photome- 
try referencing the measurements to three stars within the 
field. The corresponding lightcurve is plotted in the bot- 
tom panel of Fig. ^ The measurement errors are smaller 
than the symbol size. In addition, we include three data 
points, taken from the long-term monitoring by Schramm 
et al. (1994). Those are marked by triangles. 

2.3. Lightcurve analysis 

As evident from Fig. |l|, 0235+164 is variable in all three 
radio bands and in the optical. A major fiare around 
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Fig. 1. Intensity variations of 0235+164 in Oct. 1992 at 
1.49 GHz (A = 20 cm), 4.86 GHz (A = 6 cm), 8.44 GHz 
(A = 3.6 cm), and in the optical R-band (A = 650 nm) 
(from top to bottom). Plotted is the flux density (in Jy 
for the radio data, in magnitudes for the optical data) ver- 
sus Julian Date. For the radio lightcurves Gaussians fitted 
according to Equation (1) are included (see Section 2.3 
and Table 13). For the optical lightcurve we included three 
data points (marked by triangles) measured by Schramm 
et al. (1994). 



JD 2448905 can be identified throughout the radio fre- 
quencies, and may be tentatively connected with the op- 
tical maximum at the beginning of the observation. We 
note, however, that the exact position of the latter can- 
not be determined precisely due to the sparse sampling 
of the optical data. Therefore, we consider the connection 
between radio and optical variations as possible, but not 
definitive. 

In addition, a second flare towards the end is present in 
the 21 cm-data, possibly corresponding to the increase at 
6 cm, and the sharp peak by a factor of two in the optical. 
A corresponding feature at 3.6 cm would be expected well 
inside the observation period but is definitely not present. 
The lightcurve at 6 cm shows additional faster variations 
which have no corresponding features at the other wave- 
lengths. These faster variations (which are not shown by 
the calibrator sources and therefore are real) could for 
example be caused by scattering in the ISM as we will 
discuss later. But we note that the "global" behavior is 
very similar at all three radio wavelengths. 

We focus in this paper on the first flare (JD ^ 
2448910) which is pronounced in all three radio frequen- 
cies, and could be connected to the optical increase around 
JD 2448900. We assume that all four observed lightcurves 
are caused by the same physical event in the source. 
In order to describe this major feature, we fit a linear 
background and one Gaussian component to the radio 
lightcurves (using all data points before JD 2448910) ac- 
cording to 



S (t) ^ ao + ai ■ t + 02 • exp 



(t-asf 



(1) 



where S(t) is the measured flux density. The parameters 
and estimated errors are listed in Table nl The fits reveal 
three properties which make the radio variability quite 
unusual. First, the relative amplitude of the flare becomes 
larger with increasing wavelength. Second, the duration 
of the event (i.e., the width of the Gaussian given by the 
parameter 04) decreases with increasing wavelength. And 
third, no monotonic wavelength dependence of the time 
of the peak can be found. Including the sparse optical 
data the sequence is rather: 650 nm — > 20 cm — > 6/3.6 cm 
(the peaks at 6 cm and 3.6 cm are simultaneous within the 
errors). 

We determine the time lags between the peaks by de- 
riving Cross Correlation Functions for the radio data sets 
(again using only data before JD 2448910.0). The CCFs 
were computed using an interpolation method (e.g. White 
& Peterson 1994 and references therein). Afterwards, the 
time lags were determined by the calculation of a weighted 
mean of the CCF (i.e., the center of mass point) using all 
values > 0.5. The resulting time lags are: 



6 cm 



3.6 cm 



3.6 cm 



0.84 days, 
0.71 days, 
-0.24 days. 
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Background 




Gaussian flare 




A [cm] 


Slope 


Amplitude 


Center 


Duration 




«i [Jy/d] 


12 [Jy] 


as [JD] 


at [d] 


20 


(4.58 ±0.002) -10-^ 


0.318 ±0.01 


8904.3 ± 0.1 


1.67 ±0.1 


6 


(2.18 ±0.002)- 10^2 


0.299 ±0.01 


8905.0 ± 0.1 


2.44 ± 0.2 


3.6 


(1.99 ±0.003) -10-2 


0.243 ± 0.02 


8904.8 ± 0.3 


3.65 ±0.9 



Table 1. Fit to the radio lightcurves by a linear background and one Gaussian component 



with an error of about 0.2 days in each case. The differ- 
ences between the time lags derived from the Gaussian fits 
and the CCF are within the errors and probably due to 
the fact that the flares are not perfectly Gaussian (this ex- 
plains also that rl?^" ± t^^^ ^ Ti%lZ)- The deviation 
from the Gaussian shape is particularly obvious in the 
lightcurve at 6 cm. Nevertheless, the CCF analysis con- 
firms the result that the sequence of the flares is unusual, 
since the 20 cm peak clearly precedes the peaks in the 
other bands, while the time lag between the 6 and the 
3.6 cm-data does not appear to be significant. 

To check the significance of the time lags between the 
maxima, we carried out Monte-Carlo-Simulations for the 
Cross-Correlations between the radio frequencies. As a 
start model for the lightcurves we used the Gaussian fit pa- 
rameters with the original sampling and added Gaussian 
noise by a random process. In a second step, we allowed 
the sampling pattern to be shifted in time randomly and 
independently for every single simulation. This procedure 
confirmed that the peak at 20 cm significantly precedes 
the other two. 

2.4- Long-term variability 

In Fig. we present radio data at 4.8, 8.0 and 14.5 GHz 
obtained within the Michigan Monitoring Program (e.g. 
AUer 199S| and references therein) from January 1991 to 
November 1995. Our VLA observations (indicated by the 
arrow in the bottom panel of Fig. y) coincide with the 
peak of a large flux density outburst. 

Also in the mm- and cm-radio data, published by 
Stevens et al. ( 1994 ) a maximum at the time of the 
VLA-observations can be seen at least at 22, 37, 90 and 
150 G Hz. The optical data presented by Schramm et al. 
( |l994|) also give a clear indication for an outburst at visi- 
ble wavelengths right before our VLA observations. Three 
data points which are close to our observations are in- 
cluded in our R-band lightcurve (Fig. yj, marked as tri- 
angles. The long-term monitoring implies that our obser- 
vations took place when 0235+164 was in a very bright 
state. 

3. Discussion 

3.1. Problems 

On the basis of the collected data, and the analysis in 
the previous section, we note several properties of the ob- 
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Fig. 2. Long-term monitoring of 0235+164 since 1991 
with the UMRAO 26-meter telescope at 4.8 GHz, 8 GHz, 
and 14.5 GHz (from top to bottom). An arrow indicates 
the epoch of our VLA observations, close to the peak of 
the large flux density outburst. 

served flare which are unusual and require an explanation 
in the framework of physical models: 

— The sequence of the flares is rather unusual. The 20 cm 
maximum precedes the maxima at 3.6 cm and 6 cm. 
The first optical maximum - if connected to the radio 
events - is about four days earlier. 

— The peaks become narrower and stronger with increas- 
ing radio wavelength - a unique behavior which is not 



seen in other sources and not easily explained in any 
of the "standard" physical models. 
In case of an intrinsic origin of the variability, one can 
derive the corresponding source brightness tempera- 
ture from the duration of the event (e.g. Wagner & 
Witzel 1995). For A = 20 cm this yields Tb ~ 7- lO" K, 
far in excess of the inverse Compton limit (Kellermann 
& Pauliny-Toth|1969D. 



— Our observations show that variations are present at 
both radio and optical wavelengths, with very simi- 
lar timescales. The gaps in the optical lightcurve do 
not allow us to establish a one-to-one correspondence 
between individual events in both wavelength ranges, 
but it seems plausible that they are caused by a com- 
mon physical mechanism. This is a severe difficulty 
for models that attribute the variations to strongly 
wavelength-dependent propagation effects (free-free 
absorption and interstellar scintillation). 

In the following, we discuss various models which could 
describe the variations and take into consideration at least 
some of the peculiar properties mentioned above. 



3.2. Relativistic shocks 

Propagation of a relativistic shock-front through the jet 
is commonly accepted as one of possible causes of flux- 
density variability in AGN (e.g. Blandford & Konigl 1979, 
Marscher & Gear 1985). The time scales usually involved 
in these models are of the order of weeks to months (cor- 
responding to source sizes in the range of light weeks 
to months), and are consequently significantly longer 
than the ones observed here. Following Marscher & Gear 
( 1985 ), the characteristics of the flux density evolution 
in the case of a moving shock within the jet can be de- 
scribed as follows. Starting at high frequencies (in the 
sub-mm-regime) the outburst propagates to longer wave- 
lengths while the peak of the synchrotron spectrum shows 
a very special path in the Sm-^^m-plane. This path can be 



described by three power laws Sm oc v!^ (with different ex- 
ponents fc), distinguishing three different stages of the evo- 
lution (see also Marscher 19901 ). During the synchrotron 
or the adiabatic expansion stages, which are likely to be 
found in this wavelength range, the spectral maximum is 
expected to move from higher to lower frequencies, with 
the peak flux density being either constant or decreas- 
ing with decreasing frequency. Thus, for this "standard- 
model" , we expect that the flux density reaches its max- 
imum at higher frequencies first, and that the amplitude 
of the peak decreases towards lower frequencies. This is 
contrary to our observations. 

In contrast, the canonical behavior for a shock-in-jet 
model is seen in the long-term lightcurve (Fig. 0): The 
amplitude increases with increasing frequency, resulting 
in a strongly inverted spectrum during the outburst and 



the sequence of the peaks (determined from CCFs) follows 
the expectations: 14.5 GHz -^ 8.0 GHz -^ 4.8 GHz. 

It should be noted, that the model of Marscher & Gear 
(1985) is based on three assumptions: (i) the instanta- 
neous injection of relativistic electrons, (ii) the assump- 
tion that the variable component is optically thick at the 
beginning of the process, and (in) that the jet flow is adi- 
abatic. Therefore, this model describes a transition from 
large (r > 1) to small (r < 1) optical depths for each fre- 
quency. It is possible, however, that 0235-1-164 is initially 
optically thin at our observing wavelengths, and that the 
optical depth increases with time, e.g. due to continuous 
injection of electrons or field magnification, or through 
compression. In this case, r may reach unity - and the 
flux density its maximum - at lower frequencies earlier 



than at higher ones (e.g. Qian et al. 1996), as observed. A 
similar behavior was discussed for CTA 26 by Pacholczyk 
(1977) - although on longer time scales. In this model, 
we expect the maximum at 4.9 GHz to precede the one at 
8.4 GHz or that they are reached at the same time. The 
latter may be true within the uncertainty. However, the 
different amplitudes and the durations of the event cannot 
be explained without additional assumptions. 

Alternatively, the observed variations may be ex- 
plained with a thin sheet of relativistic electrons moving 
along magnetic field lines with a very high Lorentz factor 
(7 ~ 20-25). In this case, a slight change of the viewing 
angle (e.g. from 0° to 2-3°) may give rise to dramatic 
variations of the aberration angle and therefore of the ob- 
served synchrotron emission (Qian et al., in preparation). 
Additionally, this should cause significant changes in the 
linear polarization (strength and position angle), which 
may be studied in future observations. 



3.3. Precessing beam model 

We now investigate a scenario in which the observed effect 
is caused by the variable Doppler boosting of an emitting 
region moving along a curved three-dimensional path. If 
the observed turnover frequency of such a region falls be- 
tween 1.5 and 8.4 GHz, peaks in the lightcurves can be 
displaced relative to each other. The Doppler factor varia- 
tions required to reproduce the observed timelags may be 
caused by a perturbed relativistic beam (cf. Roland et al. 



1994, see also Gamenzind & Krockenberger 1992). The jet 

10) beam 
: 0.4. The 



is assumed to consist of an ultra-relativistic (7 : 
surrounded by a thermal outflow with speed /3 
relativistic beam precesses with period Pq and opening 
angle fio- The period of the precession may vary from 
a few seconds to hundreds of days. Roland et al. (1994) 



show that this model can explain the observed short-term 
variability of 3C 273, and also makes plausible predictions 
about the kinematics of superluminal features in parsec- 
scale jets. We use a similar approach to describe the flux 
evolution of 0235-1-164. The trajectory of an emitting com- 
ponent inside the relativistic beam is determined by colli- 
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Fig. 3. Model lightcurves at 1.4, 4.9, and i 
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mation in the magnetic field of the perturbed beam, and 
can be described by a helical path. In the coordinate sys- 
tem (x,y,z) with z-axis coinciding with the rotational axis 
of the helix, the component's position is given by 



X — r{z)cos{ujt — kz + (/jq) 
y = r(z) sin(a;t — kz + ^q) 
z = z{t), 



(2) 



where r(z) describes the amplitude of the helix, and can 
be approximated as r{z) — roz/{ao + z). For a precess- 
ing beam, ao = vq/ tanQo. The form of the function z{t) 
should be determined from the evolution of the velocity 
/3b of the relativistic component. /3b can be conveniently 
expressed as a function of z, and in the simplest case as- 
sumed to be constant. Then, under the condition of instan- 
taneous acceleration of the beam {dz/dt > 0, for z -^ 0+), 
the component trajectory is determined by 



t{z) =to + 



—dz = to 

z 



C2{Z) 



kLor-^{z) + Ci{z) 



dz , (3) 



with Ci(z) = K(z)]2 + 1, C2(z) = Ci{z) + fc2r2(z), and 
Cz{z) ^ [C2{z)(3l - cj2r2(z)Ci(z)]i/2. (This follows di- 
rectly from Equation (7) in Roland et al. (1994).) Gener- 
ally, both Lo and k can also vary. Their variations should 
be then represented with respect to z, and ijo{z) and k{z) 
be used in Equation (^. 

We describe the emission of the perturbed beam by 
a homogeneous synchrotron spectrum with spectral in- 
dex a = —0.5, and rest frame turnover frequency v[^ = 
150 MHz. The beam precession period Pq = 200 days, and 
VIq = 5?7. r{z) is described by tq = 0.1 pc and oq = Ipc. 
The corresponding lightcurves are plotted in Figure 0. 

We can see that the model is capable of reproducing 
the observed time lag between 1.4 GHz and the higher ra- 
dio frequencies. One can speculate that a more complex 
physical setting (e.g. spectral evolution of the underlying 



emission, or inhomogeneity of the emitting plasma) may 
be required for explaining the apparent discrepancy be- 
tween the modeled and observed widths of the flare. 



3.4- Free-free absorption by a foreground medium 

Here we consider the effect of free-free absorption in a 
foreground medium, either in the host of the BL Lac object 
itself, or in one of the intervening redshift systems. To 
keep the discussion simple, we neglect the cosmological 
redshift, i.e., factors (1 + z). The optical depth for free- 
free absorption of a plasma is approximately given by (see 
e.g. Lang|l974|): 



T = 8.235- lO^-'T 



-1.35 j^-2.1 



N^dl 



(4) 



where T is the electron temperature in K, i/ is measured 
in GHz, and the emission measure J N^ dl in pccm^^. 
Thus, the absorption of radiation by a foreground medium 
can be described by e"'^'^ , where c is a constant. We 
assume the following scenario. The source is moving with 
transverse speed v behind a patchy foreground medium so 
that changes in the emission measure towards the source 
produce variable absorption. To lowest order, we describe 
gaps between the clouds by 

(5) 



I 2 \2 2 
T — k X X 



{x being the axis perpendicular to the line of sight). Since 
the observed Hux density of a point source is given by 
Sohs — S e^'^ , a moving source seen through such a gap in 
the foreground medium will show peaked lightcurves with 
roughly Gaussian shape. The width of the peaks will de- 
crease with increasing wavelength. However, there are two 
major problems that need to be addressed. Firstly, in this 
model the maxima for all frequencies are reached at the 
same time, and secondly, the observed durations (i.e., the 
widths of the Gaussians fitted according to Equation (1)) 
do not follow the expected behavior 04 ex A^^. 

The time lags between the peaks of the observed 
lightcurves can be explained for an extended source by 
a slight shift of the brightness center depending on the 
frequency. 

To deal with the second problem, we assume that the 
source is not point-like, but has a circular Gaussian shape, 
with the source size proportional to the wavelength. Thus, 
the flux density is given by 

{x ~ vtf 



S{x, t) = So ■ exp 



(6) 



with ax = (7q ■ X (i.e., ctq is the source size at Im wave- 
length) . 

Assuming that the angular size of the variable region 
is much smaller than the antenna beam, the observed flux 
density is given by the integral 



Sobs(t)= / 5(x,i)-e-^(")dx 



(7) 



Evaluating the above integral, gives (since 5'obs is the 
integral of a product of two Gaussians) 

,2+2 



5'obs (i) = S' ■ exp 



-vH^ 



al + l/fc2A2 



(8) 



with a new normalization constant S'. Therefore, the 
square of the width of the Gaussian is 



2 1/2 
"4 - ^ ( f^A 



1 \ ^^^2,^1 
fc2A2 J v^ k^ W2 A2 



(9) 



and it should depend on wavelength like A- X^ + B ■ X~^. 
By adjusting the parameters A and B to fit the mea- 
sured values of al at the three observing wavelengths, 



we derive values for Uq/v and k 



We assume here 



that the transverse speed is dominated by superlumi- 
nal motion with v/c — /3app and obtain a source size 
of 0.0067 • /3app pc corresponding to an angular size of 
1-6 • /3app /^as at A = Im. We note that such a small 
source diameter even for /3app = 10 results in a bright- 
ness temperature of about 10^^ K, and therefore violates 
the inverse Compton limit. For higher velocities - as ob- 
served in this source (e.g. Chu et al. 1996) - the observed 



size can be larger. However, to reconcile our observational 
findings with the inverse Compton limit, Doppler factors 
of the order of 100 are needed. 

The second term in Equation (0) gives the size of the 
gap in the foreground medium, i.e., the distance between 
the points where t — \. Since we assumed r = k^ X^ x^, 
this distance is Ax = 2/(fc • A), which is about 2.3 • 10"'* • 
/3app pc at A = 1 m. (Note that this is true only for the case 
where the absorber is at the redshift of the BL Lac object; 
the ratio of the angular diameter distances of emitter and 
screen has to be applied as a correction factor in the case 
of an intervening absorber.) 

We still have to check whether Equation (4) gives a suf- 
ficient optical depth for reasonable choices of electron tem- 
perature and emission measure. The strongest constraints 
come from the data at 3.6 cm: to explain the observed 
amplitude of 0.24 Jy at a source flux of 5Jy, r must be 
at least 0.05 at this wavelength. For an electron temper- 
ature of 5000 K, an emission measure of 5 • lO^pccm"^ 
is needed. The thickness of the absorber cannot be much 
larger than the transverse scale derived above, which is 
0.06 pc at 3.6 cm for /3app = 10; this gives an electron den- 
sity of 10"* cm~'^. These values are within the range found 
in Galactic HII regions and planetary nebulae. 

We conclude that this model can explain the observed 
shorter duration of the flares at longer wavelengths, and 

- under the assumption of slightly different spatial loca- 
tions of the brightness center at the observed wavelengths 

- also the sequence of the peaks. It predicts that the am- 
plitude of the peaks increases more strongly with wave- 
length than observed, but it is consistent with the data 
when an underlying non- variable component is taken into 
account. However, in the possible case of a connection be- 
tween the radio and the optical variations this model fails. 



since the optical radiation would not be affected by free- 
free-scattering. 

3.5. Interstellar scattering (ISS) 

Scattering processes in the interstellar medium are well 
known to cause flux density variations at radio frequencies 
(e.g. Rickett 1990). In this section we investigate the pos- 
sibility that ISS is the cause of the variations seen in our 
observations. We will follow mainly the considerations and 
notations of Rickett et al. (1995). For a point-hke source, 
the spatial scale of flux density variations caused by RISS 
is given by (L is the path length through the medium, 
^scat the scattering angle) 



ro^x~0.25Le, 



(10) 



which is proportional to X^'^, for a Kolmogorov-type 
medium (Rickett et al. 1984), and therefore also 6'scat esc 
A2-2 (cf. Cordes et al. 1984). The spatial scale of an ex- 



tended source (assuming a Gaussian shape of ax in width) 
is then given by 



r0,\ 



O.A 



(0.5 L^a) 



(11) 



Then, the scintillation index mg_\ and the variability 
timescale Tg^\ for the extended source can be derived by 



T(9,A 



mo,x- 

re,x 

V 



' rg\ 



(12) 



where V is the velocity of the Earth (i.e., the observer) 
relative to the scattering medium, and toq.a is the (wave- 
length dependent) scintillation index of a point source. 

We assume a source diameter which is proportional to 
A as we did in the previous section, thus (J\ = ctqA, and 
use 6'scat = 6'oA2-2 (sec above). This gives 



me\ 



Te\ 



mo,\ Oq A2 



4 + 4a2A2 



^^elX^-^ + AalX- 



and 



(13) 



Therefore, it is clear that - independent of the wavelength 
dependence of mo, a - the timescales of the variations be- 
come shorter for decreasing wavelengths. This is contrary 
to our observational findings (see Table nl), implying that 
this simple model is unlikely to explain the observations. 
Additionally, interstellar scattering cannot cause variabil- 
ity in the optical regime. Hence, in this case again, a pos- 
sible connection of the optical and the radio variations 
would rule out ISS as the only cause of the observed vari- 
ability. 

However, owing to the small source diameters in- 
volved here, ISS can be present as an additional effect. 
As an example, we calculate the scintillation index and 
the timescales with the following assumptions. Following 
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Rickett ( |1986 ) the path length in the interstellar medium 
of our galaxy is L ~ 500 pc • esc |6| ~ 788 pc (the source 
galactic latitude is —40°). With ao — 1-2 mas (which cor- 
responds to Tb ~ 10^^ K), 9o ~ 60 mas and a typical 
velocity (of the observer) y = 50 km/s this yields: 




Therefore, the faster variations which are clearly seen at 
higher frequencies (especially in the 6 cm lightcurve) may 
be due to ISS. 

3.6. Gravitational Microlensing 

Another possible explanation for the origin of the ob- 
served variations is gravitational microlensing (ML) by 
stars in a foreground galaxy. ML effects have been unam- 
biguo usly o bserved in the mul tiple Q SO 2237+0305 (Irwin 
et al. 1989 , Houde & Racine 1994), and most likel y also 
in other multiply- imaged QSOs (see Wambsganss 1993, 
and references therein) . The possibility that ML can cause 
AGN variabili ty ha s long been predicted (Pacz yhsk 



1986 



Kayser et al. |1986| , Schneider & Weiss |1987[ ), but it re- 
mains unclear whether ML causes a substantial fraction 
of the observed variability in QSOs (e.g. Schneider 1993| ). 
0235+164 has a foreground galaxy {z = 0.524) situ- 
ated with in tw o arcseconds from the line of sight (Spinrad 
& Smith 1975), and an additional galaxy 0'.'5 away from 
the source (Stickel et al. 1988, see also additional com- 
ponents reported in Yanny et al. 1989). Additionally, a 
nearby absorption system was observed at A = 21cm by 
Wolfe, Davis & Briggs (1982). All three objects may host 



microlenses affecting the emission from 0235+164. Thus, 
for 0235+164 the probability for ML is expected to be 
high (Narayan & Schneider 1990| ), so that sometimes ML 
events should be present in the lightcurves. 

We will show now how ML can modulate the under- 
lying long-term lightcurve and explain faster variations of 
long-wavelength flux compared to short-wavelength radi- 
ation, even when the longer wavelength radiation comes 
from a larger source (component). Since the available data 
do not permit a detailed account of possible ML situations, 
the attention here is restricted to two simple situations: an 
isolated point-mass lens in the deflector, and a cusp singu- 
larity, formed by an ensemble of microlenses (Schneider & 
Weiss 1987, Wambsganss 1990). In fact, both cases yield 



similar predicted ML lightcurves. The scales of the source 
size and the lens mass necessary to yield a flux variation 
of the observed kind can be estimated for both cases to- 
gether. 

We assume an elliptically shaped emitting feature that 
moves relativistically in the direction roughly coinciding 
with the minor axis of the ellipse. Such a component can 



be formed by relativistic electrons which are locally ac- 
celerated by a shock front inside a super luminal jet. The 
shape of the source component and its orientation is then 
determined by the flow inside the jet. A Gaussian bright- 
ness profile is assumed, with component size ex A (see Fig. 
for details). We postulate that the emission peaks at all 
three wavelengths are displaced relative to each other, but 
that the peaks of shorter- wavelength components are situ- 
ated within the half intensity contour of longer- wavelength 
components. 

Let Pc be the apparent effective transverse velocity 
of the source component; using the redshift Zg — 0.94 
of the object, this corresponds to an angular velocity of 
Wa = 2/?/i 10~^ mas/day. If a source component moves 
along a track in the source plane, and the component size 
is much smaller than the minimum angular separation da 
from the singularity, as indicated in Fig. ^ (solid ellipse), 
then the timescale of variation is given roughly by the 
ratio da/va. On the other hand, if a strongly elongated 
source component moves so that parts of it cross the line 
of sight to the singularity (as indicated in Fig. ^, dashed 
ellipse), then the shortest possible timescale is roughly the 
ratio between the transverse angular source size a\ (the 
minor semi-axis at wavelength A) and the angular veloc- 
ity. Now assume that the former case approximates the 
3.6 cm source and the latter case approximates the 20 cm 
source. If Aia.g ^ 4 days, Aig ~ 3 days and At2o ^ 2 days 
are the variability timescales for the three wavelengths 
considered, we have 



wa Aig.e -- 8/3/1 10""* mas 



and 



r 020 ^ 4/3/i 10 mas 



(14) 



(15) 



where r < 1 is the axis ratio of the Gaussian source com- 
ponent. In order for the 20 cm source to experience ap- 
preciable variations, the closest separation of its center 
from the singularity cannot be larger than its major semi- 
axis, i.e., 020 ^ c^a, and this inequality can be satisfied for 
r < 0.5. 

Since the relative contribution of the moving compo- 
nent to the total flux of the source is unknown, we cannot 
use the observed lightcurves to determine the magnifica- 
tion of the component emission. The magnification of a 
point source at separation 6 from the point singularity is 



fJ-P 



cVJp 



(16) 



where x = 9/9o, and 6q is the angular scale induced by a 
point mass lens of mass M: 



70 = 



4GM Dds 



D.D, 



(17) 



where D^, D^, and Dds denote, respectively, the angular 
diameter distances to the lens, the source, and from the 
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Fig. 4. Geometry of the proposed microlcnsing scenario. 
Elliptical surface brightness contours are drawn for the 
3.6 cm (solid), 6 cm (dotted), and 20 cm (dashed) source 
components. The caustic point of the microlens is at the 
origin, denoted by 'X', and the source components are 
assumed to move along the solid line as indicated. The 
scales indicated on the axis correspond to the synthetic 
lightcurves shown in Fig. 0. 

lens to the source, m — M /Mq is the lens mass in units 
of the solar mass. Assuming that the lens is situated at 
z = 0.525, 



70 



= 1.87> 



mas 



(18) 



Approximating the point-source magnification by /ip ~ 
1/x (for a; <C 1), and assuming as before that the size 03. g 
is much smaller than the closest separation of the source 
from the point-like singularity, the maximum magnifica- 
tion of this source component becomes 



M3.6,i 



2Mra^''^ h-^/^ [3~ 



(19) 



Hence, a solar-mass star would yield a magnification of 
the order of 2 for the smallest source component moving 
at roughly the speed of light and in general can produce 
lightcurves similar to the observed variations. 

In Fig. 0, we plot numerically determined ML light- 
curves for a moving source with an axis ratio r = 0.4, 
minimum separation d^ — 8/3/il0~'* mas, and semi-major 
axis of the 20 cm source component of 020 = /9/il0~^ mas. 
The lens mass is m = 0Af3^h. The source sizes are cho- 
sen to be proportional to wavelength, and the brightness 
peaks of the 6 cm and 20 cm components are displaced 
relative to the peak of the 3.6 cm component by 0.4 of 
their corresponding sizes. As can be seen from the modeled 
lightcurves, the variability timescale of the 20 cm compo- 
nent is considerably shorter than that of the shorter wave- 
length components, in accordance with our analytical es- 
timates. In addition, the observed shift of the brightness 
peak at 20 cm before those at smaller wavelengths can be 



Fig. 5. Microlcnsing lightcurves obtained from the model 
discussed in the text. The three curves correspond to the 
3.6cm (solid), 6cm (dotted) and 20cm (dashed) source 
components. 

accounted for in our model by a slight tilt of the direc- 
tion of motion of the source relative to the minor axis of 
the surface brightness ellipses, in the sense of the large 
component crossing the caustic point before the closest 
approach of the 3.6 cm component to that point. Never- 
theless, we note that the small source sizes needed (in the 
range of /ias) will result in brightness temperatures of the 
order of 10^^ K, i.e., three orders of magnitude above the 
inverse Compton limit. A more detailed modeling of the 
lightcurves by a microlcnsing scenario is not warranted 
at this stage, given the large number of degrees of free- 
dom. Nevertheless, the above considerations have demon- 
strated that the basic qualitative features can be under- 
stood in the microlcnsing picture without very specific as- 
sumptions. 

4. Conclusions 

We have observed the BLLac object 0235+164 at three 
radio wavelengths and in the optical R-band and found 
rapid variations in all frequency bands. One single event 
that can be identified at all radio wavelengths shows very 
peculiar properties. The brightness peak is reached first 
at 20 cm wavelength, and afterwards at 3.6 and 6 cm. The 
amplitudes of the flares decrease from longer to shorter 
radio wavelengths, and the timescales become longer. The 
event in the radio regime might be connected to the bright 
peak in the optical lightcurve, although this connection 
remains questionable due to the sparse sampling of the 
R-band data. In the previous sections, we have discussed 
some models and to what extent they can explain the 
observed variations. 

While the conventional application of the shock-in-jet 
model has difficulties in reproducing the observations, the 
assumption of an increasing optical depth (e.g. due to con- 
tinuous injection of relativistic electrons) can cause a delay 
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of the maximum at high frequencies with respect to the 
lower frequencies, and therefore explain at least one of the 
special features. 

Variable Doppler boosting can cause simultaneous 
short-term variability in all observed wave bands. Fairly 
pronounced time lags between the different frequencies 
can be caused by turnover frequency variations in the ob- 
served spectrum of a moving source. However, broader 
peaks are expected at longer wavelengths. 

Free- free absorption and interstellar scattering are only 
capable of explaining radio variations, not variability in 
the optical regime. Therefore, if the connection between 
the optical and the radio variability is real, these mod- 
els are ruled out as the only cause for the variations. 
Furthermore, the dependence of the timescales on wave- 
length argues against an explanation of the flare by inter- 
stellar scattering. The absorption by a patchy foreground 
medium can easily describe the shape and the widths of 
the flares (in the radio) and can - if we assume different 
locations for the brightness center - also explain the time 
sequence of the brightness peaks. 

Gravitational Microlensing - in combination with a 
wavelength-dependent source size and a slight displace- 
ment of the brightness peak - provides a possible explana- 
tion for the observed variations in the radio regime. One 
would also expect fairly strong variability in the visible 
range, because of the much smaller source size. Microlens- 
ing thus appears to be a viable explanation of the observa- 
tions, which is also quite attractive because of the known 
foreground objects. 

It is quite remarkable that these attempts to explain 
the rapid radio variability in 0235-1-164 - different as 
they are - all imply that the intrinsic source size is very 
small. To reconcile the observations with the 10^^ K in- 
verse Compton limit, a Doppler factor substantially higher 
than the "canonical" value of 10 (see e.g. Ghisellini et al. 
1993, Zensus 1997) is required. Most scenarios that we 
have investigated imply V ~ 100. In this context it is 
interesting to note that circumstantial evidence for super- 
luminal motion with /?; 



app 



30 has been found in this 
source (Chu et al. 1996). The variations in 0235-1-164 are 
also among the strongest and fastest of all sources in the 
Michigan monitoring program (e.g. Hughes et al. 1992). 
This suggests that the distribution of Doppler factors in 
compact radio sources has a tail extending to 2? ~ 100, 
and that 0235-1-164 - and perhaps more generally the 
sources showing strong intraday radio variability - belong 
to this tail. The implied extremely small source size can 
allow rapid intrinsic variations, and at the same time fa- 
vor propagation effects. It is therefore plausible that the 
observed variability is caused by a superposition of both 
mechanisms. 
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